From EMC- and Cronin-effects to signals of 

quark-gluon plasma 

Wei Zhu^^, Jianhong Ruan^ and Fengyao Hou^ 

^Department of Physics, East China Normal University, Shanghai 200062, P.R. China 
^Kavh Institute for Theoretical Physics China, CAS, Beijing 100190, P.R. China 



a; 

; Abstract 

The EMC- and Cronin-effects are explained by a unitarized evolution equation, 
where the shadowing and antishadowing corrections are dynamically produced by 
gluon fusions. For this sake, an alternative form of the GLR-MQ-ZRS equation 
is derived. The resulting integrated and unintegrated gluon distributions in proton 
and nuclei are used to analyze the contributions of the initial parton distributions to 
the nuclear suppression factor in heavy ion collisions. A simulation of the fractional 
energy loss is extracted from the RHIC and LHC data, where the contributions of 
the nuclear shadowing and antishadowing effects are considered. We find a rapid 
crossover from week energy loss to strong energy loss at a universal critical energy 
of gluon jet Ec ~ lOGeV. 
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1 Introduction 

One of the important findings at RHIC and LHC is tliat liigli transverse momentum kf 
hadron production in central heavy ion coUisions is suppressed compared to p+p coUisions 
[1, 2]. This suppression can be attributed to energy loss of high-Zcj partons that traverse 
the hot and dense medium (i.e., quark-gluon plasma QGP) formed in these collisions. An 
important goal of the study of heavy ion collisions is therefore to determine the properties 
QGP from the measured fractional energy loss after deducting nuclear effects on the initial 
parton distributions. 

The parton densities in a bound nucleon differ from that in a free nucleon. One of 
such examples is that the ratio of nuclear and deuterium's structure functions is smaller 
or larger than unity at Bjorken variable x < 0.1 or 0.1 < x < 0.3. These two facts are 
called as the nuclear shadowing and antishadowing in the EMC effect [3]. The nuclear 
shadowing and antishadowing effects origin from the gluon fusion (recombination) between 
two different nucleons in a nucleus, which changes the distributions of gluon and quarks, 
but does not change their total momentum [4] . In consequence, the lost gluon momentum 
in the shadowing range should be compensated in terms of new gluons with larger x, 
which forms the antishadowing effect. 

An other example is the Cronin effect: the ratio of particle yields in d + A (scaled by 
the number of collisions)/^ + p, exceeds unity in an intermediate transverse momentum 
range (Cronin enhancement) and is bellow unity at smaller transverse momentum (anti- 
Cronin suppression). This effect was first seen at lower fixed target energies [5] and has 
been confirmed in ^/s — 200GeV d + Au collisions at the BNL Relativistic Heavy Ion 
Collider (RHIC) [6]. 

The Cronin effect is more complicated than the EMC effect, the former mixes the 
shadowing- antishadowing corrections at initial state and the medium modifications at 



final state. The later is an important information for understanding the properties of 
dense and hot matter formed in high-energy heavy-ion collisions. Therefore, the nuclear 
shadowing and antishadowing effects, which have appeared in the EMC effect, should be 
extracted from the Cronin effect, then the remaining result exposes the medium properties. 

The saturation models are broadly used to study the Cronin effect. The satura- 
tion is a limiting behavior of the shadowed gluon distribution in the Jalilian-Marian- 
lancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation [7], where the uninte- 
grated gluon distribution is absolutely fiat in fc^-space at kt < Qs, Qs is the saturation 
scale. An elemental QCD process, which arises nonlinear corrections in the JIMWLK 
equation is also the gluon fusion gg g. As we have pointed out that the antishadowing 
effect always coexists with the shadowing effect in any gluon fusion processes due to a 
general restriction of momentum conservation [8]. However, such antishadowing effect is 
completely neglected in the original saturation models. The Cronin enhancement in these 
models, (i) is additionally explained as multiple scattering [9] using the Glauber-Mueller 
model [10], or the McLerran-Venugopalan model [11]; (ii) is produced by special initial 
gluon distributions of proton and nucleus [12]. A following question is: what are the 
nuclear antishadowing contributions to the Cronin effect? 

A global Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) analysis of nuclear par- 
ton distribution functions (for example, the ESP09-set [13]) was proposed, where the data 
from deep inelastic scattering (DIS), Drell-Yan dilepton production, and inclusive high-Zcj 
hadron production data measured at RHIC are used. They found that a strong gluon 
antishadowing effect is necessary to support the RHIC data. However, the shadowing 
and antishadowing effects in the DGLAP analysis are phenomenologically assumed in the 
initial conditions, since the DGLAP equation [14] does not contain the nonlinear cor- 
rections of gluon fusion. Due to the lack of the experimental data about nuclear gluon 



distribution, the above mentioned status are undetermined. For example, a similar global 
DGLAP analysis shows that the available data are not enough to fix these comphcated 
input distributions and even appearing or not the antishadowing effect are uncertain 
[15]. Besides, the DGLAP equation in the coUinear factorization scheme evolves the inte- 
grated parton distributions, the behavior of the unintegrated gluon distributions, which 
contain information of the transverse momentum distribution is completely undefined in 
this method. Therefore, the ESP09-set of nuclear parton distributions can not predict 
the RHIC data at lower kt, where the contributions from intrinsic transverse momentum 
become important [13]. 

The modification of the gluon recombination to the standard DGLAP evolution equa- 
tion was first proposed by Gribov-Levin-Ryskin and Mueller-Qiu (the GLR-MQ equation) 
in [16,17], and it is naturally regarded as the QCD dynamics of the nuclear shadowing 
since the same gluon fusion exists both in proton and in nucleus, and they differ only 
by the strength of the nonhnear terms [18]. However, the GLR-MQ equation does not 
predict the nuclear antishadowing effect since the momentum conservation is violated in 
this equation. This defect of the GLR-MQ equation is corrected by a modified equation 
(the GLR-MQ-ZRS equation) by Zhu, Ruan and Shen in [19,20], where the corrections 
of the gluon fusion to the DGLAP equation lead to the shadowing and antishadowing 
effects. Unfortunately, the integral solutions in the present GLR-like equations need the 
initial distributions on a boundary line {x, Qq) at a fixed Qq, and they still contain the 
unknown input with nuclear shadowing and antishadowing effects at small x. 

Motivation of this work is trying to improve the above questionable methods. We shall 
study the nuclear shadowing and antishadowing effects in the EMC- and Cronin-effects, 
which are dynamically arisen from the gluon recombination. Then we use the resulting 
nuclear gluon distributions to predict the contributions of the initial parton distributions 



to the nuclear suppression factor in heavy ion coUisions and to extract fractional energy 
loss from the RHIC and LHC data. For this sake, an alternative form of the GLR-MQ-ZRS 
equation at the double-leading-logarithmic-approximation (DLLA) is derived in Sec. 2. 
The evolution of this equation is along small x-direction. The nonlinear corrections to the 
input distributions may neglected if the value of the starting point Xq is large enough. The 
shadowing and antishadowing effects are naturally grown up in the evolution at x < Xq. 
Thus, we avoid the unperturbative nuclear modifications in the input distributions and 
simplify the initial conditions. We use the existing data about the EMC- and Cronin- 
effects to fix a few of free parameters in the solutions. Then the resulting integrated 
and unintegrated gluon distributions in proton and nuclei are used to analyze the nuclear 
suppression factor in heavy ion collisions. 

Our main conclusions are: (i) we support the stronger shadowing-antishadowing effects 
in the heavy nucleus both in the unintegrated and integrated gluon distributions due to 
a strong ^-dependence of the nonlinear corrections; (ii) the anti-Cronin suppression and 
Cronin enhancement mainly origin from the same gluon recombination mechanism as in 
the nuclear shadowing and antishadowing effects of the EMC effect; (iii) fractional energy 
loss is rapid crossover from week energy loss to strong energy loss at a universal critical 
energy of gluon jet ~ lOGeV. 

The organizations of this work are as follows. We shall derive the GLR-MQ-ZRS equa- 
tion in a new form in Sec.2. Using this equation we study the shadowing and antishadow- 
ing effects in the EMC effect in Sec. 3. These resulting unintegrated gluon distributions 
in proton and nuclei are used to expose the nuclear shadowing and antishadowing contri- 
butions to the Cronin effect in Sec.4. The nuclear shadowing and antishadowing effects 
to the nuclear suppression factor are predicted and a simulations of fractional energy loss 
is extracted from the RHIC and LHC data in Sec. 5. 



2 A new form of the GLR-MQ-ZRS equation 

In history, the DGLAP evolution is derived by using the renormahzation group method 
for the integrated distributions. The resulting equation evolves with factorization scale 
II. In this section, we try to write the DGLAP equation with the nonlinear modifications 
beginning from the unintegrated distributions. Then, we get an alternative form of these 
equations, which evolves with the Bjorken variable x. 

We begin from a deeply inelastic scattering process, where the unintegrated gluon 
distribution is measured. In the /c^-factorization scheme, the cross section is decomposed 
into 



da{probe*P kX) 
( k"^ X \ 

= /(xi, k\^ (g) K, — , ® da(probe*ki k) 
\kit xi J 

= Af{x, kl) da{probe*ki k), (1) 

which contains the evolution kernel /C, the unintegrated gluon distribution function / 
and the pro6e*-parton cross section d(T{probe*ki — > k). For simplicity, we take the fixed 
QCD coupling at the leading order (LO) approximation in this work. According to the 
scale-invariant parton picture of the renormalization group [21], we regard A/(a;, /c^) as 
the increment of the distribution /(xi, k\^ when it evolves from (xi, k\^ to (x, k^). Thus, 
the connection between two functions f{xi,kl^) and f{x,kf) via Eq. (1) is 

f{x,e,) = f{x^,kl,) + Af{x,e,) 

In the case of the LO DGLAP evolution, we adopt a physical (axial) gauge, which sums 
over only the transverse gluon polarizations, so that the ladder-type diagrams dominate 



the evolution, the unintegrated distributions satisfy the normahzation relation 



G{x,i,') ^ xg{x,i,') = J ^xf{x,k^) ^ J ^F{x,k',), (3) 

t,min t t,min t 

where the possible non-logarithmic tail for kt > n are beyond NLO accuracy. These distri- 
butions correspond to the density of partons in the proton with longitudinal momentum 
fraction x, integrated over the parton transverse momentum up to kt = ji. 
Prom Eqs. (2) and (3), we have 




where the last step is held when the /c^-strong ordered. Usually DGLAP evolution is 
written in terms of the virtuality k"^ rather than k^ , but at LO level this is the same since 
the difference is a NLO effect. Therefore, we take 

G(x,e,)^G(x,i,'). (6) 

Thus, in 

G(x,i,')^G(x,,i,l) + AG(x,i,'), (7) 

we write 
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at the leading logarithmic approximation, in which the unregularized splitting kernels 
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We add the contributions of the nonlinear recombination terms in Eq. (8) according 
to Refs. [19, 20] (See appendix A), 
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G'(^)(x,/x^)=it:GG^(x,//^), (12) 

where = 1/ (ttT?^) is a correlation coefficient with the dimension [L""^] , it! is the effective 
correlation length of two recombination gluons. One can easily get the GLR-MQ-ZRS 
equation at DLL approximation 



dG{x,i^^) 
din fi'^ 

TT Jx Xi 7lfJ,^R^ — 1 Jx Xi 

^n^m^N^-lJx/2 ^^(^^'^)- (13) 
It is interesting to compare this small-x version of the GLR-MQ-ZRS equation with 
the GLR-MQ equation, which is [17] 
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where 



G^'\x,^^')^^^G\x,^^% (15) 

is assumed. 

Comparing with the GLR-MQ equation (14), there are several features in the GLR- 
MQ-ZRS equation (13): (i) the momentum conservation of partons is restored in Eq. 
(13); (ii) because of the shadowing and antishadowing effects in Eq. (13) have different 
kinematic regions, the net effect depends not only on the local value of the gluon distri- 
bution at the observed point, but also on the shape of the gluon distribution when the 
Bjorken variable goes from x to x/2. In consequence, the shadowing effect in the evolu- 
tion process will be obviously weakened by the antishadowing effect if the distribution is 
steeper. Therefore, the antishadowing effect can not be neglected in the pre-saturation 
range. 

According to the definition Eq. (3), one can roughly estimate the unintegrated gluon 
distribution using 



^dG{x,id^) 
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(16) 



However, Eq. (16) cannot remain true as x increases, since the negative virtual DGLAP 
term may exceed the real emission DGLAP contribution and it would give negative values 
of F. In fact, due to strong kt ordering in DGLAP evolution, the transverse momentum 
of the final parton is obtained, to leading-order accuracy, just at the final step of the evo- 
lution. Thus, the /ct-dependent distribution can be calculated directly from the DGLAP 
equation keeping only the contribution which corresponds to a single real emission, while 
all the virtual contributions from a scale equal to kt up to the final scale ji of the hard 
subprocess are resummed into a Sudakov factor T. The factor T describes the probability 
that during the evolution there are no parton emissions. However, at the small x range, 
the virtual contributions to the gluon distribution in the DGLAP kernel can be neglected. 
Besides, we have indicated that the contributions of the virtual processes in the nonlin- 
ear kernels of the GLR-MQ-ZRS equation are canceled [19], therefore, the Sudakov form 
factors are same in nucleon and nucleus [22] and they are really canceled in the following 
ratio form. Thus, we use 
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and obtain 



^F{x,kl) = II- 



lt,rmn LZ -L 

Now we re-derive the GLR-MQ-ZRS equation, but which evolves with the longitudinal 
momentum. We differentiate 



F{x, kl) = F{x,, kl) + AF(a;, k^), (19) 
with respect to x. Note that 
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Generally, the resummation solution is hard to obtain from this equation. However, at 
the LL{l/x) approximation the second term on the right-hand side of Eq. (20) vanishes. 
In this case, the resummation becomes simple, i.e., we have 
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, (0.15 <x< 0.3), 
(21) 

This is a new form of the GRL-MQ-ZRS equation. The negative and positive nonhnear 
terms correspond to the shadowing and antishadowing effects in the gluon recombination. 
Note that the shadowing and antishadowing coexist in the region x < 0.15, while there is 
only the antishadowing in 0.15 < x < 0.3. 

The explanation of the kinematic regions of Eq. (21) is follows. The evolution kernel 
of the GLR-MQ-ZRS equation (11a) is derived at LL((5^) approximation as same as the 
derivation of the DGLAP equation (9a) and they are valid in a whole x range. However 
the DLLA form of the GLR-MQ-ZRS equation (lib) works at the small x range, for say 
X < xi, where a;i = C(10~^) according to ln(l/a;) ln(A;|.//x^) ~ We take xi = 0.15 

in this work. 

Now we apply Eq. (21) in the nuclear target. Shadowing and antishadowing effects 
are thought to arise from a nonlinear mechanism which occurs when gluons are sufficiently 
dense to interact themselves. The strength of the gluon recombination is proportional to 
the gluon density on the transverse area. The gluons with smaller x exceed the longitudi- 
nal size of nucleon in a nucleus. Thus, the strength of the nonlinear recombination terms 
in Eq. (21) should be scaled by A^^^ in a nucleus. On the other hand, although the softer 
gluons of different nucleons with extra small kt maybe correlated on the transverse area 
because the integrations on kt can go down to a very small value in Eq. (21), we neglect 
these corrections due to F{x, /c^ ^ 0) — )■ 0. In this simple model, Eq. (21) in the nucleus 
becomes 
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The input gluons distribute on the boundary hue (2x0, /ct) at fixed 2xo. We take a 
larger value of xq as a starting point of the evolution, where the gluon fusion just begins 
and the nonlinear corrections to the input gluon distributions can be neglected. On the 
other hand, the contributions of the Fermi motion and nuclear binding effects to the 
nuclear parton distributions at small x are small. In this nucleus is simply the 

incoherent sum of the constituent nucleons at 2xo and we have 



F{2xo,kl) = FAi2xo,k^), (23) 

where the nuclear parton distributions have been normalized by derived A. 

The DGLAP kernel in Eq. (21) resums the double leading asl'n.{l/x)ln{k^ / fi^) con- 
tributions. As we know that the BFKL evolution [23] which resums the leading ln(l/x) 
contributions and should completely replace the DGLAP equation at vary small x range, 
for say at a; < 10"^ according to agln^l/x) ~ In this regime, Eq. (2) is rewritten 

as 
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In consequence, we have the hnear BFKL equation (See Appendix B) 
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and nonhnear BFKL equation with the modifications of gluon fusion 
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To make our progress we need to know how the BFKL formahsm hnks with the con- 
ventional DGLAP dynamics. A usual framework which gives a unified treatment through 
out the {x, kt) kinematic region has been provided by Catani, Ciafaloni, Fiorani and 
Marchesini (the CCFM equation [24]). This equation is based on the coherent radiation 
of gluons, which leads to an angular ordering of the gluon emissions along the chain. In 
the leading ln(l/x) approximation the CCFM equation reduces to the BFKL equation, 
whereas at moderate x the angular ordering becomes an ordering in the gluon transverse 
momenta and the CCFM equation becomes equivalent to standard DGLAP evolution. 

Unfortunately, the CCFM schema contains the unknown shadowing and antishad- 
owing information in its complicate input distributions. In this work we use following 
physical picture. As we know that the BFKL equation can be derived in a dipole picture 
(see Appendix B). At lager x region (for say, x > 0.1), parton densities in nucleon are 
dilute and the probe is interacted with a single parton (Fig.l a). In this case, the DGLAP 
dynamics are dominated. On the other hand, at very small x region {x < 10~^), the 
correlations among the initial partons in a nucleon becomes important and the dipole 
configuration dominates the initial state (Fig.lc), and the BFKL dynamics replace the 
DGLAP equation. Note that although x < 10~^ is a dominated regime by the BFKL 
evolution according to Q;sln(l/a:;) ~ 0{1), we can not exclude the BFKL dynamics par- 
ticipate the part of evolution at a more larger x region since asln{l/x)ln{k^/ii^) < 1 
allows X > 10~^. Therefore, a natural connection between two evolution dynamics is that 
the BFKL dynamics replace asymptotically the DGLAP dynamics from x = 0.1 to 10~^ 
through the mixing of the single parton and dipole configuration as shown in Fig. lb. In 
the concretely, we take 
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We emphasize that if we use the original form of the GLR-MQ-ZRS equation (13) to 
replace Eq. (21), the solution of the Eq. (31) becomes very difficult since existing two 
different evolution ways. 

All parameters in the solutions of Eqs. (21), (22) and (29) will be determined by the 
EMC effect. Most data about the EMC effect are measured by the structure functions. 
Thus, we should calculate the quark distributions at small x. We assume that the sea 
quark distributions at the small x range are dominated by the gluon distribution, via the 
DGLAP splitting g qq. Thus, the deep inelastic structure function at small x reads 



[27] 




(32) 



where Pqg{xi) is the DGLAP splitting function. 



3 The EMC effect 

We choose xq = 0.15 as the stating point of the evolution in Eqs. (21) and (22), where 
the nonlinear gluon recombination begins work. We find that following input is satisfy, 
i.e., 

F{2xo, kl) = 2^)exp{-{log{kl)f). (33) 

The computations of Eqs. (21), (22) and (29) need pre-know the value of F[xi/2, k^) 
at the step x — Xi. For this sake, we proposed following program in [26] 

fXi 2\ TP f^i ,2\ , ^linear {f,kt) - Fshadowing (f^^t) 

V 2" *y ^Shadowing [^^'^tj^ iA - A + 1 ' ^ ^ 

where Fshadowingi^i/'^, kj) (or Fiineari^i/'^-i kj)) indicates that the evolution from Xi to Xi/2 
is controlled by Eqs. (21), (22) and (29) but without the antishadowing contributions (or 
is controlled by the linear equation). The parameter < A < 1 implies the different 
velocities approaching to the dynamics of Eqs. (21), (22) and (29). 

At first, we use the well known F2p{x, (5^)-data [25] of a free proton to determine the 
parameters in the computations. Then we predict the distributions in nuclei. The dashed 
curve in Fig. 2 is our fitting result using the input (33), where we take the parameters 
R = 2.4 GeV-^, k^ = 0.01G'ey^ a, = 0.3 and A = 0.02. Note that the contributions 
from the valence quarks to F2 at x > 0.1 are necessary and they can be parameterized by 
the difference between the dashed curve and experimental solid curve in Fig. 2. 

Figure 3 shows our predictions of Eqs. (21), (22) and (29) for the Ca/C, C/Li, Ca/D 
and Cu/D compared with the EMC and NMC results [28,29]. The agreement is accept- 
able. 

It is different from the scheme of the DGLAP evolution equation, we can predict the 
nuclear effects for the unintegrated gluon distribution. The results are presented by the 



ratios of the nuclear unintegrated gluon distributions in Figs. 4 and 5. Although deep 
inelastic scattering experiments do not provide the direct test of these effects, hadron- 
nucleus scattering at RHIC relates to the nuclear unintegrated gluon distributions and we 
will use them in next section. 

The ratios Gca{x,Q'^)/GD{x,Q'^) for gluon distributions at = 2 and 10 GeV"^ 
using Eqs. (21), (22) and (29) are given in Fig.6. The Q^-dependence of the gluon ratio 
is predicted in the region 10~^ < x < 10"^ in our model. The logarithmic slope b in 
Ga/ Ga' = a + blnQ^ is positive. However, the corresponding slope in the ratio of the 
structure functions -F2Ca/-^2D is negative at small x (see Fig.7). For example, the predicted 
(5^-slope for calcium at a; ~ 4 x 10~^, b ~ —0.046, and at x ~ 10^^, b ~ —0.03, the results 
are compatible with the measured data in [30] . 

The data of Gsn{x, Q^)/Gc{x, Q^) are measured from inelastic J/ip production by the 
NM Collaboration in Ref. [31], which determine a stronger nuclear antishadowing effect 
but with a larger uncertainty. Our prediction is presented in Fig.8. 

Comparing with the ESP09 set [13], we predict a more stronger antishadowing effect in 
the gluon (integrated and unintegrated) distributions of heavy nucleus. One of the reasons 
is that the observed antishadowing effect in the structure functions origins dynamically 
from the gluon fusions in our model, while in the DGLAP analysis it is partly from the 
input distributions of the valence quarks [13]. The second reason is that the A-dependent 
nonlinear terms enhance the effect of the gluon fusion in the heavy nuclei. 



4 The Cronin effect 



The Cronin effect is described by the nuclear modification factor RdA, which is defined 
as the ratio of the number of particles produced in a, d + A coUision over the number of 
particles produced in ap + p collision scaled by the number of collisions 
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kt and r) are the transverse momentum and the pseudo-rapidity of the observed hadron, 
respectively. Ncoii is the number of collisions in d + A scattering. In Eq. (35), 
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where z — kt/kt^g] Dp{z) and Da{z) are the fragmentation functions of gluon jets in 
proton and nucleus, where the factorized scale-dependence of the fragmentation functions 
are neglected; The rapidity y of the produced gluon in the center-of-mass frame of p -|- p 
collisions is defined by 



2^1/2 = •exp(±7/); 



(38) 



The relation between the rapidity y and pseudorapidity 77 for massive particles is 



where me// is the typical invariant mass of the gluon jet. 

We assume that the hadrons in the central region are produced from the hadronization 
of the gluons in gg ^ g mechanism. According to Ref. [16] we have 



dap-i,{h-t. q) 
(Pktdr) 

ANr r dz a. 



r/=0 



fdz as 2n f \ 

J ' dqlgtpl{x,{kt,g-qt,gf)ipl{x,qlg) 
ANc [dz as 2r> i ,\ 



kt,g=Jkt/z 



(^(x, kl)G^{x, kl) + G^{x, kDip^ix, kl) 



kt,g=Jkt/z ' 



(40) 



dad-A{kt,rj) 



d^ktdr] 
ANr r dz a 



r)=0 



Vc [dz as 2n ( \ 

rk^ 

J dqlg^lix, {kt,g - qt,gf)ip^{x, qlg) 
ANr r dz a 



kt,g=Jkt/z 



-J 



Z Kt,g 



roAz) 



'<Pl{^. kL)G^{x, kl) + G'^ix, kl)ip^{x, kl 



kt,g=Jkt/z ' 



(41) 



■•gy-^^'^gJ^ \-^^'H,g/ ' ^ '%gJ^g '%gJ 

where ipg{x,ql) = F{x,qlg)/qfg and we neglect the /c^-dependence in the fragmentation 
functions. 

At first step, we neglect the interactions at final state, i.e., in Eqs. (40) and (41) 



Dp{z)^Da{z)^S{z-1). 



(42) 



We indicate this nuclear modification factor as wliicli is drawn in Fig. 9. Ac- 
cording to Eq. (38) at y = and ^/s — 200Gey, the antishadowing effect on the gluon 
jet should distribute in a broad range SGeV < kt < GOGeV, which corresponds to the 
antishadowing range 0.02 < x < 0.3 in Fig. 4. 

In next step, we consider the corrections from the fragmentation functions but ne- 
glecting the difference between proton and nucleus. We take [32] 



where D{z) ^ at z ^ since the coherence effects in QCD at small z. Our results are 
plotted by the sohd curve in Fig. 10. The data are taken from the BRAHMS results in 
[6]. One can find that the fragmentation functions shift cross point between the Cronin 
and anti-Cronin effects toward to small kt, since the position of the peak value of the 
fragmentation functions always localizes at small z. We find that the nuclear shadowing 
and antishadowing effects at the initial state dominate the Cronin effect, although a small 
nuclear modification to the fragmentation functions, i.e., Da{z) ^ Dp{z) is necessary. 

The parton energy loss caused by medium-induced multiple gluon emission in various 
nuclear conditions is a hot topic, since the data in the RHIC Au-\-Au collisions show a new 
hot matter and it is maybe a strongly coupled Quark-Gluon Plasma (sQGP). The presence 
of a dense medium influences the space-time development of the partonic shower of a jet. 
An energetic jet parton propagates through the medium, a part of its energy transfers 
to the thermal partons, and it is called the parton energy loss. After this jet propagates 
a longer distance in an expanding de-confinement system, most gluons carrying the lost 
energy escape from the jet cone, and they are un-measured. Thus, for an example, the 
modified fragmentation function in an effective model can be written as [33] 




(43) 



Da{z) 



1 




z 



(44) 



1-e 



where E is the initial energy of a gluon jet and e = AE/E is fractional energy loss. 

We consider that a similar energy loss mechanism also exits in the RHIC d + Au 
coUisions but with a smaller value of e. The dashed curve in Fig. 10 is the resulting 
nuclear suppression factor with e = 0.1. From the above mentioned results we find that 
the anti-Cronin suppression and Cronin enhancement origin from the nuclear shadowing 
and antishadowing effects in the initial state as in the EMC effect. 

Using the above mentioned parameters we predict the Cronin effect in d + Pb collisions 
at the CERN Large Hadron CoUider (LHC). The results in Fig. 11 show that the anti- 
Cronin suppression dominates a broad kt range, while the Cronin enhancement appears 
at large kt range. This character is different from the RHIC data and it comes from the 
kinematical x range. The value of x in LHC is much small, where the nuclear shadowing 
dominates the gluon distribution. While the value of x in RHIC can reach the nuclear 
antishadowing range. 



5 The signals of QGP 

One of the important findings at RHIC and LHC is tliat liigli transverse momentum 
hadron production in central heavy ion coUisions is suppressed compared to (properly 
scaled) p+p collisions [1, 2], which is defined as the nuclear suppression factor Raa- 
This suppression can be attributed to energy loss of high-Zcj partons that traverse the 
hot and dense medium formed in these collisions. However, the extraction of the energy 
loss from Raa needs to know the nuclear effects in the initial parton distributions. Our 
obtained integrated and unintegrated gluon distributions in proton and nuclei provide 
such information. 

Some of the most important results obtained at RHIC are related with the high-Zc^ 
spectrum in heavy ion Au + Au collisions. In this aspect, a precise determination of 
the nuclear effects in the initial state of the collision is fundamental. In what follows we 
present the estimation for the ratio Raa 



dNA-A(kt,ri) 
cPkTdr] 



^^coll (pktdrt 



(45) 

??=o 



using our parameters in the explanation of the EMC- and Cronin-effects. At first step, 
we assume that e in Eq. (44) and show the nuclear shadowing and antishadowing 
effects in the nuclear suppression factor. The result is dashed curve in Fig. 12. One 
can find that a big difference between dashed curve and data at kt > 5GeV, which is 
commonly interpreted in terms of a strong energy loss of the energetic partons when 
traversing a dense medium. For example, we take e — OA and the result is the pointed 
curve in Fig. 12. Obviously, a true form of fractional energy loss e is crossover from the 
smaller value of e the larger one when increasing energy of the gluon jet. It is interest 
that we take 



e — 



a if E <Ec 
b if E> Er 



where a — 0.2 and b — 0.4 (i.e., Fig. 13a); similar to Eq. (41) we have 
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ifEc{l-b)<kt<Ec{l-a); 

i/ /ct > £;c(i-a), 



kt,g=Jkt/z 



kt,g=Jkt/z 



kt,g=Jkt/z 

(47) 



we find that the resulting (solid curve) in Fig. 12 perfect to fit the data at IGeV < kt < 
lOGeV. 

Recently, the nuclear modification factor in central Pb+Pb collisions at = 2.76TeV 
is published by the ALICE Collaboration at the LHC [2]. The data indicate that Raa 
reaches a minimum of about 0.14 a,t kt — 6 — IGeV , which is smaller than at RHIC, 
however, it steeply rises over the asymptotic value of RHIC at larger kt- Therefore, it 
is unclear whether a more denser mater is formed at LHC. Obviously, a quantitative 
determination of the energy loss excluding gluon shadowing and antishadowing effects is 
necessary. For this sake, similar to the above mentioned approach in the RHIC data, we 
take a = and b — 0.58 (Fig. 13b) in Eq. (46) and the result is presented the solid curve 



in Fig. 14. The dashed and pointed curves correspond to e = and 0.58, respectively. 
Although the data show that the energy loss e will decrease with the jet energy E » 
10 GeV, Eq. (46) is a good approximation describing the nuclear suppression factor at 
IGeV < kt < lOGeV. Thus, we find a rapid crossover from week energy loss to strong 
energy loss at a universal critical energy of gluon jet ~ lOGeV. Using same parameters 
we predict the nuclear suppression factor in central Pb + Pb collisions at y/s = 5.5 and 
7.7TeV in Figs. 15 and 16. 

In summary, the EMC- and Cronin-effects are explained by a unitarized evolution 
equation, where the shadowing and antishadowing corrections are dynamically produced 
by gluon fusions. For this sake, an alternative form of the GLR-MQ-ZRS equation is 
derived. The resulting integrated and unintegrated gluon distributions in proton and 
nuclei are used to analyze the contributions of the initial parton distributions to the 
nuclear suppression factor in heavy ion collisions. A simulation of the fractional energy 
loss is extracted from the RHIC and LHC data, where the contributions of the nuclear 
shadowing and antishadowing effects are considered. We find a rapid crossover from week 
energy loss to strong energy loss at a universal critical energy of gluon jet Ec ~ lOGeV. 
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Appendix A. GLR-MQ-ZRS equation: The modifications of tlie gluon recombination 
to tfie DGLAP evolution in tfie GLR-MQ-ZRS equation lias following form [19,20]. 



dG{x,Q^) 
(ilnQ2 



Ja:/2 \Xi J 



for gluon distribution and 



dxS{x, Q^) 



Ja::/2 Xi Xi \Xi / 

for sea quark distributions, where P^^ are the evolution kernels of the linear DGLAP 
equation and the recombination functions 



dXi GG^GG, LL{Q'^) 
/^GLR-MQ-ZRS 

Xi 

—f^ 5 -dxi {A - 3 

64 XiX 



dXi ^GG^SS, LL{Q^) 
'^GLR-MQ-ZRS 

a'^J2xi-x)^(18x'i -21xix + 14:x) , 

li- — h -dxi {A - 4) 
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Appendix B. A derivation of the BFKL equation in the dipole picture: 

In order to certify the physical picture in Fig. 1, we give a derivation of the BFKL 
equation in the dipole picture. We consider the DIS process in Fig. 17, where the dashed 
lines are the time-ordering lines in the TOPT. These processes imply a scattered gluon 
before its radiation is omitted from two correlating gluons. We call such a correlating 
gluon cluster as a dipole, which phenomenologically describes the correlation among initial 
partons. 

The evolution kernel in a QCD evolution equation is a part of more complicated 
scattering diagram. In general, the correlations among the initial partons make these 
partons off mass-shell. One of the most important hypotheses in the derivation of the 
DGLAP equation is that all correlations among initial partons are negligible during the 
interaction. Therefore, the interaction of a virtual probe with the nucleon can be factorized 
as the nonperturbative parton distribution and hard probe-parton matrix in the coUinear 
factorization scheme. At the higher density range of partons, these correlations among 
the initial partons can no longer be neglected. In this case, the transverse momenta of 
the initial partons are non-zero and these partons are generally off mass-shell, therefore, 
the /cT-factorization scheme is necessary. 

In this work we use the semi-classical Weizsacker- Williams {W — W) approximation 
[34] to realize the /cj-factorization scheme. The reason is that the W — W approximation 
allows us to extract the evolution kernels and to keep all initial and final partons of the 
evolution kernels on their mass-shell. 

For convenience, we use k to indicate the transverse momentum kt in this section. 
We consider that the observed wave function 'if{x2,k) is evolved from the initial wave 
functions and ^(xi,^^) via the QCD interactions, i.e.. 



where two perturbative amplitudes corresponding to Fig. 17 are 



(B-l) 



BFKLi = M ^ \ ^ , ^ —M^, {B-2) 



and 



' ' M2. (5-3) 



V ^Pa + ^Pb Ek + Ei^ — Ep^ 

The momenta of the partons are parameterized as 



Pa - {x,P + ^1^, ^ + L, xiP), {B 4) 



k^{x2P + ^,k,X2P), {B-5) 

2X2 P 



la - {{X^ - X2)P + ^, (a^l - a;2)P), {B - 6) 

2(Xi - X2)P 



and 



The matrixes of the local QCD interactions are 



Ml = igf "[gaoiPa + k)^ + gp^{-k + + gja{-la - Pa)/:i]ea{Pa)ep{k)e^{la), {B - 9) 



where the polarization vectors are 



<^) = (0,e,-|^), {B-12) 

and 

where e is the transverse polarization of the gluon in = (eo,e, £3) = (0,e, 0), since the 
sum includes only physical transverse gluon states in the TOPT form [19]. 

Taking the LL(l/x) approximation, i.e., assuming that X2 <^ a^i, one can get two 
similar amplitudes 



Abfkli = tgr'''2M^, {B - 14) 

V X2 k 

and 



Abfkl2 = igr''2.l^^^. {B - 15) 

However, these two amplitudes really occupy different transverse configurations. This is 
a reason why the dipole model of the BFKL equation is derived by using the transverse 
coordinator-space. However, we shall show that the momentum representation still can 
be used to distinguish the differences between Eqs. (B-14) and (B-15). 



The two parton correlation function is generally defined as 



= /(^'^c,k6), (5-16) 

where and fc^^ are conjugate to the impact parameter and transverse scale of a cold spot. 
Equation (B-16) implies the probability finding a gluon, which carries the longitudinal 
momentum fraction a; of a nucleon and locals inside a dipole characterizing by and k^f,. 

In this work we derive the evolution equations in the impact parameter-independent 
case. This approximation implies that the evolution dynamics of the partons are domi- 
nated by the internal structure of dipole. Thus, the evolution kernel is irrelevant to k,. 
and we shall use 

f{x,k^) = J ^f(x,k,,k^), {B - 17) 

which has the following TOPT-structure 



f{x,kab) 
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ab 



2Ef 



M, 
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n 2 



n 2 



2E, 



ab 



y (27r)32^x' 



{B - 18) 



.Ep — Eab — Ex- 

Notice that all transverse momenta in Eqs. (B-1)-(B-15) are indicated relative to 
the mass-center of the nucleon target. However according to Eq. (B-17), the evolution 
variable is the relative momentum therefore, it is suitable to rewrite all momenta to 
relative to in Eq. (B-2) and to in Eq. (B-3), respectively. Thus, we replace the 
transverse momenta as follows: 



and 



in Eq. (2.6) since 



and 



La ~^ hah bob — P „ ~ kt — haOi 



hab — kaO + 



{B - 19) 



{B - 20) 



and 



Lb ~^ kba ~ koa — 2ib~ k — kbo, 

in Eq. (B-3). In consequence, we have 



{B - 21) 



and 



{B - 22) 



[xl 




i_ kob 


/ X2 




k^ 



{B - 23) 

where we identify two e in Eq. (B-23) since the measurements on (2:2,^0) {x2,k^ij) 
are really the same event. 

Equation (B-1) with Eqs. (B-22) and (B-23) provide such a picture: a parent dipole 
with the longitudinal momentum fraction xi and transverse momentum radiates a 
gluon, which has the longitudinal momentum fraction X2 and the transverse momentum 

^0 (oi" kob)- 



We take the square of the total amphtude, one can get 



E, 



ab 



2Ef 
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SE Epj-Q^^ 
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{2ny2Ek' 



= f{xi,kib) ® —JCbfkl {kab,kaa, o^s) ® da{q*^„f^^k{x2,kxo) ^ k'{x2,k!)) 

= A[^{x2,k^oW{x2,kao) + '^{x2,kaoW{x2,koh) + 



Af(x2,k^o) (g) da(q^„,,^k(x2,kao) k'(x2,^)), 



{B - 24) 



where on the last step, because the probe only picks up the contributions from ^'(x2, ^ao)**(^2, 
we regard Af{x2,k^o) as the increment of the distribution f{xi,kab) when it evolves from 
(3^1,^6) to {x2,kao). Therefore we have 



A/(a;2,fc«o) 



dkab dXi Xi 



k^^ Jx2 Xi X2 



K.BFKL {tab, koQ, as) f{xi,kib), 



{B - 25) 



or 



AF(a;2,^o) = ^X2f{x2,k^o) 

= [ I —^BFKL{kab^ka()^Ois)F{xi,kab)- 



{B - 26) 



Using definition 



-^(2^2,^0) ^ F(xi,k^i,) + AF(x2,kao), 



{B - 27) 



we write 



dF{x,k^o) 



dx 

= j d^b^BFKLikab-, kaO, Ols)F{x, kab) , {B - 28) 

According to Eq. (B-22), the evolution kernel reads as 



X2 Xi ZXi [ZTT) 

otgNr k^,h dxi ,9 , , „ , 

- I2 .2 ,2 -^d'k.b- {B-29) 



Finally Eq. (B-28) becomes 



^ dF{x,k^Q) 
dx 

= ^ / d-'k^-S^Fix^k.,). {B - 30) 

This is the real part of the BFKL equation. 

Using the TOPT-cutting rules [19], one can prove that the virtual diagrams in Fig. 
17 contribute the similar evolution kernel as the real kernel but differ by a factor — 1/2 x 
(1/2 + 1/2). Here we call the cut diagram with a naive partonic amplitude without any 
QCD-corrections as the virtual diagram. The negative sign arises from the changes of 
time order in the energy denominators. The factor (1/2 + 1/2) is due to the fact that the 
probe "sees" only the square root of the parton distribution accepting the contributions 
of the partonic processes in a virtual diagram, and the other factor 1/2 origins from the 
symmetry of the pure gluon process. Therefore, the evolution equation corresponding to 
Fig. 17 is 



^ dF{x,k^b) 
dx 



= / ^^^0 ,2,, \ (i^ - 31) 

Since we calculate the contributions to AF{x, k^o), we should make the replacement 6 -H- 
in Eq. (B-30). Combining the real and virtual parts of the evolution equation, we have 



^ dF{x,k^o) 
dx 



(fk 



^ J 2 7,2 F{^ykab) ~ 7,2 7,2 -^(^)^ao) 



{B - 32) 



27r2 

According to Eq. (B-18) , the distribution f{x, k) in the TOPT-form contains a singular 
factor 1/k^, which arises from the off energy-shell effect in the square of the energy 
denominator. In order to ensure the safety of taking the W — W approximation, we 
remove this factor to the evolution kernel and use the following new definition of the 
unintegrated gluon distribution 



F{x,k) = ^F{x,k), (S-33) 
k 

where is a unity vector on the transverse momentum space. Thus, Eq. (B-32) becomes 



dF{x, k^o) 



O^s^c f ,2 7, k.aO 



dx 

9^2 / d%b-rl^ U - F{x, fc„o)] , {B - 34) 

which consists with a usual form of the BFKL equation. 

The correlations among the initial gluons can be neglected in the dilute parton system. 
In this case the contributions of the interference diagrams Figs. 15c and 15d disappear. 
Thus, the kernel (B-29) reduces to the sphtting functions in the DGLAP equation at the 
small X limit. 



xidxi 2n OsNcdxid^ 

X2 Xi 77 X2 k 



= ^DGLAP 




{B - 35) 



The two initial gluons have the same transverse momentum and we always can take 
it to zero and use the coUinear factorization to separate the gluon distribution. The 
corresponding DGLAP equation reads 



where the scaling restriction S{x2 — xb) is included and Eq. (3) is used. 

The above derivations reveal a relation between two evolution equations from a new 
view point: the BFKL equation dominates the evolutions at the small x regime, where 
the initial partons appear in the dipole picture. On the other hand, the parton densities 
become dilute at the larger x regime. In this case, the dipole decomposes to the isolate 
partons and the BFKL kernel reduces to the DGLAP kernel due to the interferant part 
of the BFKL kernel disappears. 





{B - 36) 



Appendix C: 

We compare the solutions of the following DGLAP equation 



dx 

lt,mirt J-C 

with that of the BFKL equation (28) using the same input (33), to illustrate that the hnear 
terms of Eq. (21) is a reasonable form of the DGLAP equation at the DLL approximation. 
The sohd and dashed curves in Figs. 18 and 19 show the unintegrated gluon distributions 
in the DGLAP (C-1) and BFKL (28) equations, respectively. One can find that the 
emissive gluons are random walk in /c^-space in the BFKL evolution, while they are ordered 
in kf in the DGLAP evolution. In consequence, the distribution steeply increase at small 
X in the BFKL equation more than that in the DGLAP equation. 
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{a)x > 0.1 (6)10-3 < a; < 0.1 {c)x < lO'^ 



1 The illustrating kinematic regions of the DGLAP- and BFKL equations. 




Fig. 2 The fit of the computed F2p(x, — 10 GeV"^^ in proton by the evolution 
equations (21), (29) and (32) using the input Eq. (33) (dashed curve). The contributions 
of the valence quarks are parameterized by the differences between solid and dashed 
curves. The data are taken from Ref. [27]. 




Fig. 3 Predictions of the evolution equations (21), (22) and (29) compared with the 
EMC ratio of the structure functions for various nuclei. The data are taken from [28,29]. 
All curves are for = 10 GeV'^ 




Fig. 4 Predictions of the evolution equations (21), (22) and (29) for the ratio of the 
unintegrated gluon distributions in Ca/D with different values of x and given kt. 




5 Similar to Fig. 4 but with different values of kt and given x. 




Fig. 6 x-dependence of the ratio for the integrated gluon distributions inCa/D in the 
evolution equations (21), (22) and (29). 




7 Similar to Fig. 6 but for the the ratio of the structure functions. 




Fig. 8 Predictions for the ratio of the gluon distributions in Sn/C and the data are 
taken from Ref. [31]. 




Fig. 9 Predicted nuclear modification factor R^^^u gluon jet in central d + Au 
collisions at \/s = 200GeV. 




Fig. 10 Nuclear modification factor RdAu oi charged particles in central d+Au collisions 
at ^/s = 200GeV, where fractional energy loss e = (solid curve) and 0.1 (dashed curve). 
The data are taken from Ref . [1] . 




Fig. 11 Predicted nuclear modification factor RdPb of charged particles in central 
d + Pb collisions at (a) ^/s = b.bTeV and (b) 8.8TeV. 




Fig. 12 Estimated nuclear suppression factor R^a in central Au + Au collisions at 
^/s = 200GeV: solid curve using Eq. 46 with a = 0.2, b = 0.4, (see Fig. 13a); dashed 
curve using e = 0, and pointed curve using e = 0.4. The data are taken from Ref . [1] . 
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Fig. 13 Two possible fractional energy losses, which correspond to (a) central collisions 
at ^ = 200^6 V for Au + Au and (b) at = 2.76TeV for Pb + Pb, respectively. 
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Fig. 14 Similar to Fig. 12, but for central Pb + Pb collisions ect ^/s ^ 2.76TeV, where 
solid curve using Eq. 46 with a = 0, b = 0.58, (see Fig. 13b); dashed curve using e = 0, 
and pointed curve using e = 0.58. The data are taken from Ref. [2]. 




Fig. 15 Predicted nuclear suppression factor Raa in central Pb + Ph collisions at 
= 5.5TeV, all parameters are same as Fig. 14. 
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Fig. 16 Predicted nuclear suppression factor Raa in central Ph + Ph collisions at 
= 7.7Tey, all parameters are same as Fig. 14. 
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Fig. 17 The TOPT-diagrams, which lead to BFKL equation. The dashed Unes are 
the time ordered hnes in the TOPT. 




Fig. 18 The A;|-dependence of the unintegrated gluon distribution with different values 
X in the hnear part of Eq. (21) (sohd curves) and BFKL equation (dashed curves). 




19 Similar to Fig. 18 but for the x-dependence with different values of kt 



